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Abstract. The thermal regime in planetary regoliths involves three processes: 
propagation of visible radiation, propagation of thermal radiation, and thermal 
conduction. The equations of radiative transfer and heat conduction are formulated for 
particulate media composed of anisotropically scattering particles. Although the equations 
are time dependent, only steady state problems are considered in this paper. Using the 
two-stream approximation, solutions are obtained for two cases: a layer of powder heated 
from below and an infinitely thick regolith illuminated by visible radiation. Radiative 
conductivity, subsurface temperature gradients, and the solid state greenhouse effect all 
appear intrinsically in the solutions without ad hoc additions. Although the equations are 
nonlinear, approximate analytic solutions that are accurate to a few percent are obtained. 

Analytic expressions are given for the temperature distribution, the optical and thermal 
radiance distributions, the hemispherical albedo, the hemispherical emissivity, and the 
directional emissivity. Additional applications of the new model to three problems of 
interest in planetary regoliths are presented by Hapke [this issue]. 


1. Introduction 

The temperature distribution in the upper layers of the re- 
golith of a planet is important for two reasons: this region is the 
one sampled by remote measurements of the thermal radiation 
emitted by a planet, and the temperature of these layers con- 
stitutes the boundary condition for thermal models of the 
deeper interior. However, few of the published thermal models 
for this region treat radiation adequately. The classical analysis 
by Wesselink [1948] neglected both optical and thermal radia- 
tion except as a boundary condition, as did Spencer [1990] in 
his important discussion of the effects of surface roughness. 
Usually, when thermal radiation is included its only effect is 
assumed to be a term proportional to T 3 , where T is the 
absolute temperature [e.g., Henderson and Jakosky, 1994] 
(however, notable exceptions are papers by Conel [1969] and 
B. Henderson and B. Jakosky (Near-surface thermal gradients 
and mid-IR spectra: A new model including scattering and 
application to real data, submitted to Journal of Geophysical 
Research, 1 [1995]). This assumption is a good approximation 
for the deep interior, where thermodynamic equilibrium pre- 
vails, but is incorrect near the surface, where the radiation 
escapes to space and is not in equilibrium with the medium. 
Although the solid state greenhouse model of Brown and Mat- 
son [1987] allows visible sunlight to penetrate into the medium, 
it does not include the effects of scattering; hence the discus- 
sions in their paper are incomplete. This paper uses the heat 
equation and the equation of radiative transfer to describe the 
flow of energy within a particulate medium. The effects of both 
visible and thermal radiation on the temperature distribution 
are included, in addition to ordinary solid state conductivity. 
Thus radiative conductivity and the solid state greenhouse 
effect appear intrinsically in the model. 

In section 2 the basic equations that govern the thermal and 
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radiative conditions are introduced. For maximum generality, 
the quantities that appear in these equations are then reduced 
to dimensionless form. The time-independent equations are 
solved in section 3 for a particulate medium df finite thickness 
whose lower boundary is held at a fixed, temperature, and in 
section 4 for a semi-infinite medium in equilibrium with inci- 
dent sunlight. The distribution of temperature and radiance 
with depth is calculated for several representative cases. As 
with my previous radiative models, quantitative analytical so- 
lutions are derived which, although approximate, retain the 
essential physics and are sufficiently accurate for most plane- 
tary remote sensing applications. 

Only time-independent problems are treated in this paper. 
Although the differential equations for the general time- 
dependent case are given, their solutions will be considered in 
a separate paper (B. Hapke, manuscript in preparation, 1996). 
Using this model, three problems of interest in planetary re- 
mote sensing are treated by Hapke [this issue]. 

2. Basic Equations 

2.1. Assumptions and Notation 

The system that will be considered in this paper is a semi- 
infinite, horizontally stratified, particulate medium in a vac- 
uum, in equilibrium either with sunlight incident on its upper 
surface or with a lower boundary held at a fixed temperature. 
The particles of the medium are assumed to be irregular in 
shape, randomly positioned and oriented, with their light scat- 
tering and thermal properties independent of depth. An exact, 
quantitative description of the propagation of radiation 
through such a medium is beyond present capabilities, even 
with the aid of modern high-speed computers. Hence a num- 
ber of simplifying assumptions are necessary. The most impor- 
tant assumption is that the radiation fields in the medium can 
be described by the equation of radiative transfer. In principle, 
this equation applies only to radiation propagating through a 
continuous medium whose elements do not interact coher- 
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Figure 1. Schematic diagram of the geometry of the model. 


ently. However, Hapke [1981, 1993a] has argued that in a 
random particulate medium most coherent effects are negligi- 
ble and has shown that solutions of the equation provide rea- 
sonable descriptions of the bidirectional reflectances of plan- 
etary regoliths in the visible. Goguen [1996] has compared 
exact and approximate solutions of the radiative transfer equa- 
tion with measurements on a medium of microspheres whose 
scattering properties are known and found reasonable agree- 
ment. Moersch and Christensen [1995] also found that models 
based on the radiative transfer equation predicted the spectral 
emissivity of powders reasonably well. Hence this assumption 
has empirical support. 

The notation and nomenclature used in this paper are sim- 
ilar to those of Hapke [1993a], but with necessary additions and 
modifications. For convenience and reference, quantities are 
listed in the notation section at the end of the paper. A sub- 
script A denotes that a quantity is wavelength-dependent, while 
a subscript v or T denotes that the quantity is averaged over 
the visible or thermal infrared spectral region, respectively; an 
asterisk denotes a reduced quantity. 

The reduced quantities are defined as follows: 

E\{z) - £E a (z), E* v (z) = &JLz), 

E* r (z) = &E T (z). 

q = ((T 0 hT)T 4 r /E* T kT r = (T () TfhTEU 

(see below). 

R = drVdr* T = E*JE* T = E£-JE t & 

(see below). 

T r is the reference temperature that characterizes the par- 
ticular system. 

T*( t* t , t*) = T(t t , t)/T r . 

The diffusion time is t D = pCIE^k 
t* = t/t D . 


w*K = (1 - 0>a/£L < = (1 - 

W * = (1 — 

<P*,(t* v , t*) = <p v ( t„, t)/J v , t*) = < p T (r, t)/(o- 0 Tt/'ir). 

Acpt/rt, t*) = A «p„(T„, t)/J v , 

A<p*(r*, t*) = A cpj( t, t)/(cr 0 T 4 / tt). 
yl = (1 - w\Y n = Ta /£ a , y\ = (1 - w*„) 1/2 = yJC vt 
y* T = (1 - w* T ) U2 = y T /Cr- 
p a = C\p-> P* = Clp> Pt = CtP- 

p tlpo, p*)v = ClPo> Pot ~ £tPo- 
dr* = i\dr t dT*,= £ 2 v dt v , dr%= fidr-r. 

If visible irradiance is incident, choose T r as the blackbody 
radiative equilibrium temperature, T r — (i r/ ? ,/er 0 ) 1/4 , so that 
q = JJE* T kT r = < T 0 Jl /4 hrE* T k ; if/„ is the solar irradiance, 
then q = a {) ^ 3/4 / 7 rE* T kD 3/2 , where is the solar constant 
and D is the distance to the Sun in AU. If the particle prop- 
erties do not change with depth, so that E*JE* T is constant, 
then R = at all optical depths. It is usually more 

convenient to choose r* T as the depth variable, but t* t and r* 
may be readily interconverted using R. 

2.2. Equations 

The geometry is shown in Figure 1. The wavelength- 
dependent equation of radiative transfer for light that has been 
scattered one or more times within the medium, and including 
both external and thermal emission sources, is 

di K (z, a, t) 

cos'd -£ A (z)/ A (z, a, t ) 

dz 

+ j 1 - / A (z, a', t)G K (z, a', a) dtv 

4tt I 

J 4 IT 

+ g a (z. n,„ n) 

47 T 

+ — U x [T(z, t)l (1) 

7T 

where 7 x (z, a, t ) is the spectral radiance (power per unit area 
per unit steradian per unit wavelength interval at elevation z 
traveling in direction a at time t and wavelength A), d is the 
angle between a and the z axis, E x (z) is the volume extinction 
coefficient at elevation z and wavelength A, G A (z, a', a) is 
the volume coefficient at location z and wavelength A for scat- 
tering radiance from direction 7r — a' into direction a, J K is 
the spectral irradiance from the source, F(t) describes the time 
dependence of the source irradiance, p Q (t) is the cosine of the 
angle of incidence, K x (z) is the spectral volume absorption 
coefficient at location z and wavelength A (by Kirchhoffs law, 
equal to the emissivity), and U X (T) is the Planck function for 
the power emitted per unit area per unit wavelength interval 
from a blackbody at temperature T, 

U x (T(z, /)] = [2irh 0 cJ\ 5 ][e hMT{z ^ - l] -1 , (2) 


exp 


Po(t) 


£ a (z') dz' 
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where h 0 is Planck’s constant, c 0 is the speed of light in vac- 
uum, k 0 is Boltzmann’s constant, and T(z, t) is the absolute 
temperature. 

In addition to the radiative transfer equation, the tempera- 
ture must also satisfy the heat equation 

dT(z,t) d \ dT(z,t) 

PC dt dz dZ 


transfer equation (1), which can be solved to give the radiance 
at any wavelength. 

The equations will be further simplified by introducing wave- 
length-averaged parameters as follows. Integrating the radia- 
tive transfer equation over the visible portion of the spectrum 
gives 

d /visible I \{Z, Cl, t) dk 

cos-u r- 



* 

1 f * 

+ 

K x (z)J x F(t) exp 

A 

-,o(0 E ' izUz 

J 2 


K x (z) I K (z, Cl, t ) dd dk 




U x [T(z, /)] dn dk, 


( 3 ) 


where p is the bulk density, C is the specific heat per unit mass, 
k is the thermal conductivity, and K x (z) is the spectral volume 
absorption coefficient. The parameter k is the intrinsic, solid- 
state thermal conductivity of the medium only and does not 
include radiative heat transfer, which is specifically accounted 
for by the thermal radiance. 

The term on the left-hand side of (3) describes the rate of 
change of heat content per unit volume in a layer of thickness 
dz at depth z. The first term on the right-hand side of (3) is the 
difference between the heat conducted into a layer at depth z 
and the heat conducted out. The second term on the right 
describes the source of heat due to unscattered light from the 
source that has penetrated to z and been absorbed there. The 
third term is the scattered radiance absorbed at depth z, while 
the fourth term describes the loss of heat by thermal radiation 
from the particles in the element. 

The heat equation is coupled to the radiative transfer equa- 
tion only through integrals over wavelength, rather than the 
spectrally resolved radiance. In a medium illuminated by a 
source of visible light, such as the Sun, the radiance is appre- 
ciable in only two spectral regions, the visible (including the 
near-UV, visible, and near-IR) and the thermal IR, and the 
temperature distribution is determined by the integrals of the 
spectral radiance over those two regions. In many problems of 
interest the two wavelength intervals are sufficiently separated 
that the area of overlap contributes little to the integrals. This 
assumption will be made here. Note that this assumption does 
not preclude treating regions of the spectrum where both vis- 
ible and thermal radiation are important in the spectral radi- 
ative transfer equation [e.q., Hapke, 1993a, b], only that such 
regions contribute little to the integrals over wavelength. 

The system of coupled partial differential equations will be 
solved by the following procedure. The equation of radiative 
transfer will be converted into two differential equations by 
integrating over the visible and thermal wavelength regions 
separately. Since the resulting equation governing the visible 
radiation does not contain the temperature, it may be solved to 
give the visible radiance as a function of z. This solution will be 
inserted into the heat equation. The coupled heat equation 
and thermal radiative transfer equation will then be solved 
simultaneously to give the temperature and thermal radiance 
distributions. Once the temperature is known, it can be in- 
serted into the source function U X (T) in the spectral radiative 


E x (z)I x (z, Cl, t ) dk + 


Since the visible radiance consists of scattered light from the 
source, it has a spectral distribution similar to that of the 
source J x but modified by departures from a gray body caused 
by absorption bands in the medium. In order to make the 
equations tractable, it will be assumed that to a sufficient ap- 
proximation, the various parameters in the visible radiative 
transfer equation may be replaced by their averages over the 
visible portion of the spectrum weighted by the spectrum of the 
source. Thus, for example, the first term on the right-hand side 
of the integrated equation can be written 


E x (z)I x (z, a, t ) dk = E v (z) / A (z, a, t) dk, 

,| e *' visible 

where E v (z) is the average visible extinction coefficient, 
E v (z) = f E x {z)I x (z, a, t) dk I f /*(z, a, t) dk 


E x {z)J x dk 


J K dk. 

J A 


Similarly, the IR radiance at any location in the medium has 
a spectral distribution similar to the Planck function of a black- 
body at the local temperature, modified by any absorption 
bands in the medium at thermal IR wavelengths. Again, in 
order to make the problem tractable, it will be assumed that 
the parameters in the IR radiative transfer equation may, to 
sufficient approximation, be replaced by their averages over the 
IR portion of the spectrum weighted by the Planck function of 
a body at a reference temperature T r that is characteristic of 
the system. For instance, in a medium illuminated by visible 
radiation, T r might be the radiative equilibrium temperature, 
or if no visible source is present, as in a laboratory experiment, 
T r might be the temperature of the source supplying heat to 
the medium. Thus, when the radiative transfer equation is 
integrated over the thermal IR, the extinction coefficient is 
replaced by its average thermal value: 


E t (z) = 


£ a (z)/ a (z, n, o dk 


I \{z, Cl, t ) dk 


E x (z)U x (T r ) dk 


U x (T r ) dk. 


Other parameters that appear in the wavelength-integrated 
visible and thermal-IR equations of radiative transfer will be 
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replaced by their similarly averaged values. This procedure will 
not introduce serious errors if the reflectance spectrum of the 
medium is relatively bland in the visible and thermal IR near 
the maxima of J K and U K (T r ). However, if broad spectral 
features of large contrast are present there, the spectrum of 
I K (z, ft, t) could be appreciably different from J A or U x (T r ), 
which could cause wavelength-averaged parameters to be un- 
representative of their true values. 

In this paper it will be assumed that wavelength averaging 
each parameter, as described above, is valid. Let / „(z, ft, t ) 
and I T (z, ft, t) be the wavelength-integrated visible and ther- 
mal-IR radiances, respectively. Let a subscript v on any other 
quantity denote the average value of that quantity weighted by 
the spectrum of the source of visible radiation integrated over 
wavelength. Let a subscript T on any other quantity denote the 
average value weighted by the Planck function of a blackbody 
at some temperature that characterizes the system integrated 
over wavelength. 

Then the equation of radiative transfer for visible radiation 
becomes 


cos# 


dl v (z, ft, t ) 
dZ 


= —E v (z)i v (z, a, t) 


+ 4 ~ | f v (z, ft', t)G v (z, ft', a) dn r 


+ J v F(t) ^ G„(z, fto, ft) 


+ 7—1 I T (z, ft', t)G T (z, ft', ft) d£l' 
4 TT 


(To 


+ K t (z) ^ T\z, t). 

7 r 

Let t t = /” E r (z') dz' be the optical depth for IR radiation. 
Putting G r (z, ft', ft) = S T (z)p T (z, ft', ft), w T (z) = 
S t (z)/E t (z), % t (z) = K t (z)/E t (z); dividing by E T (z); 
and making t t the spatial variable, this equation can be written 
in the form 

dI T (T T , ft, t) 

cosd — ar T — = n> l) ~nr 


ft * t)p T {j T , g') dCl 


+ % T (r T ) ^ T\ r T , t), 

TT 


( 6 ) 


where S T (z) is the thermal scattering coefficient, p T (z, ft', 
ft) is the thermal phase function, w T (z) is the thermal single 
scattering albedo, and % T (z) is the volume emissivity. 
Similarly, the heat equation is 


P C 


dT(z, t) 
~dt 


dz 


dT(z, t ) 


dz 


+ K t (z) j It(z, ft, t) dCl 

477 


exp 


Pu (t) 


E v (z') dz' . 


K t (z ) I — T\z, t) da 

TT 

477 


Let t„ = J” E v {z') dz' be the optical depth for visible light. 
Putting G v (z, ft', ft) = S v (z)p v (z, ft', ft), w v (z) = 
S v (z)/E v (z), dividing by E v (z), and making r v the spatial 
variable, this equation can be written in the following form: 

dl v ( T„, ft, t) 

COS# = -/„( Tv, ft, t) 

a T v 

W„(t„) 


477 


I v (t v , ft', t)p v (r v , g') da' 


+ J v F(t) p v { r v , g 0 )e 


(4) 


where g' is the angle between it - ft' and ft, g 0 is the angle 
between ft () and ft, 5„(z:) is the visible scattering coefficient, 
p v (z, ft', ft) is the visible phase function, and w v (z) is the 
visible single scattering albedo. 

Since 


+ K v {z)J v F(t ) exp 


lio(0 


+ K v (z) | I v (z, ft, t) da, 

477 


which can be written 


dT(z, t ) d 


dt 


dz 


dT(z, t ) 


dz 


E v (z') dz' 


+ K T (z) j It(z, ft, /) da 

4tt 


O’o 


K t (z) — T 4 (z, t)] da + K v {z)JF{t)e~ r "' 

TT 

4,7 


+ K v (z) I v (z, ft, 0 da. 


(7) 


U X {T) d\ = o- 0 T 4 , 


(5) 


where cr 0 is the Stefan-Boltzmann constant, the wavelength- 
integrated equation of radiative transfer for thermal radiation 
is 

dI T (z, ft, t) 

cos# = ~E t {z)I t (z, ft, t) 

dz 


Equations (4), (6), and (7) are the fundamental equations to 
be solved. It will be assumed that all of the quantities in these 
equations, except the radiances and temperature, are indepen- 
dent of z and T. This system of equations will then be solved 
using the two-stream method. In this approximation, the radi- 
ative transfer equation is integrated with respect to solid angle 
separately over the upward going (0 < # < tt/2) and downward 
going (tt/2 < # < 7 r) hemispheres, and each quantity in each 
integral is replaced by its average value in that hemisphere. See 
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Hapke [1993a] for details. Anisotropic scattering is taken into 
account by the hemispherical asymmetry parameter 0. This 
quantity is defined such that a fraction (1 + 0)/2 of the light 
incidence on a particle from any direction within an upward 
going or downward going hemisphere is treated as scattered 
uniformly forward into the same hemisphere into which the 
light was traveling, and a fraction (1 — 0)/2 treated as scattered 
backward uniformly into the opposite hemisphere. 

Let /„,( t„, t) and I v2 ( t v , t) be the hemispherically aver- 
aged upward going and downward going visible radiances, re- 
spectively, 


I vi (An 0 2^ 


I" ir/2 


/„( t„ n, o <m ; (8a) 


J — f 


U T„ 0 = 2S n> 0 <ift - 


l=ir/2 


(8b) 


The result of hemispherical averaging of (4) for the visible 
radiation consists of two coupled equations: 


1 a/„,( r„, 0 

2 3t„ 


w, 


1-^(1 + 0») 




+ y (1 - MUr„, 0 + VW 4^(1 - W 


and 

i a/^(T„, o r , vv t , 


+ 


dr,, 


1 - 2“ (1 + 0„) 




W 

+ y (1 - () + /,/«) - 4 - (1 + 




(9b) 


i 0 _ W) 


W « l ,2 ■ P* 


Similarly, by letting 


- 2tt 


C ir/2 




(12b) 


/ r (T r , 12, t) dft, (13a) 


/ 77 ( t 7-, /) — 2^ j ^r( T r» 12, /) <ifl, (13b) 

<Pt(t t , t) = [I n ( T t, t ) + / n (Tr, 0]/2, (14a) 

A<p r (T r , t) = 0] ~ ItA'Tti 0]/2, (14b) 

and assuming that the particles of the medium radiate isotro- 
pically, the IR radiative transfer equation can be put into the 
form 


1 3A<P7-(T7-, t) 

2 d Tj 


— (1 — W t )<Pt(Tt, t ) + ^7- — T 4 (t 7-, t), 

(15a) 


- \ d(pr ^ Tr, Jl = _(i - p T w T )k<pTir T , 0, (l 5b ) 

2 dr T 

Similarly, letting y T = (1 - tv r )' /2 and = (1 — 0 r H> 7 -) ,/2 
and realizing that % T = 7 r, equations (15a) and (15b) can be 
put into the form 


1 d(p T {T T , t) 

A<p r (r r , t) - 2 ^2 dTr > 


(16a) 


and 


where the factor of ± !/2 on the left-hand side of (9a) and (9b) 
is the average value of cos# in each hemisphere. Let 

() = «) + Ut„ t)V 2; (10a) 

A<p„(t., () = <) - /„ 2 (r„, 1)1/2; (10b) 

c p v (t v , t) is the spherically averaged radiance in the medium at 
optical depth t„ and time t. By alternatively adding and sub- 
tracting (9a) and (9b), they can be put into the form 


1 aWrnO = frf ^ t) _ ^2 5 r(Trj () . (16b) 

4 OTy- 77 

Finally, the heat equation can be written in the form 


dT(z, t ) d 

pC dt = Jz 


8T(z, t ) 


dz 


+ E 1 {z)y%4 / n 


<Pt(z, t) ~ — T *( z ’ 0 


1 aA<p„(T», f) 

2 dr v 


-(1 - W„)(p v ( t„, t) +J v F{t) 


w 

Z1L e -T*f wit) 

4tt 

(11a) 


and 


1 dy ; ,(T„, t ) 

2 dT v 


0„w y )A (p v {(p v , t ) 




4ir 


0„e 


“T,/ (!))(/) 


(lib) 


Let y v = (1 - wj ,/2 and - (1 - 0 y w t ,) 1/2 . Then (lib) 
can be solved for A<p„, 


Acp^T,,, 0 


1 dcpjr,,, f) 

2^ ar. 


Wv Pv _ 


/ ” /? ( / ) 4tt Cl 


T „/n(|(() 


(12a) 


Differentiating (12a) with respect to r i; and substituting the 
result into (11a) gives 


+ E v (z)y 2 x [J v F(t)e- TjluM + 4t r<p„(z, 01* (I 7 ) 

In order to solve (12), (16), and (17), six boundary condi- 
tions must be specified. Three are that the temperature and 
radiances must be finite everywhere, including z -> 00 if the 
medium is infinitely thick. If the medium is of finite thickness 
then conditions on the lower boundary must be given (see 
section 3). Two more result from requiring that there be no 
sources of visible or IR radiance above the t = 0 level, except 
for the collimated irradiance J v , which is included as a source 
term in the equations, and is not a boundary condition. This 
requires 


o' 

II 

o' 

(18a) 

<Pv( 0, 0 - A<p„(0, /), 

(18b) 


and 
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MO, 0 = 0, (19a) 

or 

<p r (0, t ) = A<p r (0, t ). (19b) 

The remaining boundary condition on T comes from the re- 
quirement that the only flux of heat crossing the r = 0 bound- 
ary is carried by thermal radiation and none by conduction. 
However, the conducted heat flux, kdT/dz, must be continu- 
ous across any horizontal plane, including the upper surface of 
the medium. Since k = 0 above the surface, but not below it, 
this requires that 

dT 

(0, 0 = 0 (20) 

as the surface is approached from below. 

(These boundary conditions at the surface of the medium 
apply only to bodies without atmospheres. If the medium is in 
an atmosphere, then at the surface the heat flux conducted 
through the medium must be continuous with the flux con- 
ducted through the air above the surface. Also, the downward 
radiances I v2 and I T2 at the surface must be equal to the visible 
and IR radiances scattered and radiated downward by the 
atmosphere.) 

An additional requirement on the solutions arises from 
physical considerations. Deep within a medium, far from any 
boundaries, the thermal radiation must be in equilibrium with 
the temperature of the material there. This requires that 

<Pt(z, t) ^»C7 0 r 4 (z, t)/TT. (21) 


2.3. Equations in Reduced Form 

In order to make these equations as general and simple as 
possible, it is desirable to reduce them to nondimensional form 
by scaling all quantities to reference quantities that character- 
ize their magnitudes for each specific case. This has already 
been done in the radiation equation for the distance variable 
by replacing it with the optical depths t v or t t , which are 
distance in units of the extinction lengths, L Ev = E~ A or 
L et — E e 1 . A natural quantity to use to scale the visible flux 
is the incident irradiance The temperature may be scaled to 
some reference temperature T r that characterizes the system, 
and the thermal flux may be scaled to the thermal radiance 
a 0 Tf/iT inside a cavity at that temperature. If the source of 
visible radiation is sunlight, it is convenient to choose T r as the 
blackbody radiative equilibrium temperature given by J v = 
ct 0 T^/tt. Although the period of rotation P of a body might be 
used to scale the time variable, the equations become simpler 
if the diffusion time t D = pC/£ E E E k is used instead. 

In addition, if the particles of the medium are not isotropic 
scatterers, the equations may be further simplified by using 
similarity relations. These are reduced single scattering albe- 
dos and optical depths that convert the equations for anisotro- 
pic scatterers into expressions whose mathematical forms are 
as close as possible to those for isotropic scatterers. 

In the remainder of the paper an asterisk will denote a 
reduced quantity. These are defined in section 2.1. Substituting 
them into the equations formulated in 2.2 gives the reduced 
equations. For visible radiation, (12a), (12b), and (18) become 


i aVtCT*., f) , 2 


dT, 


= y :>:(<, /*) 


Fit) 


w* v 1 + pJ2p, a v* 


4tt 


1-/3, 




(22a) 


i a<p*»( t*„ c) m-: e. 






(22b) 


with the boundary condition 


(22c) 

For thermal radiation, (16a), (16b), and (19) become 
1 f*) = n _ r ,4 (T . r> ( 23a) 

1 (*) 


dT r 2 


A<P*( r' n t*) = 2 
with the boundary condition 


dT- 


1 d<p*(0, t*) 
<pr(0, t*) = 2 


(23b) 


(23c) 


Arbitrarily choosing t* t as the reduced spatial variable, the 
reduced heat transfer equation (17) becomes 

dT*(T* T , t*) d 2 T*(T* T , t*) 


dt S 


dT * T 2 


+ 4irqy T * 2 [<p* T (T* T , f*) - T* 4 (r* r , f*)] 

+ Rqy*, 2 [F(t)e~ M,,) + 4ir t*)], 
with the boundary condition 

dT*{ 0, t*)/dT* T = 0. 


(24a) 


(24b) 


Inserting (23a) into (24a) gives an alternate form of the re- 
duced heat equation: 


dT*( t t , t*) d 2 T*(T%, t*) 


dF 


dT% 2 


+ Trq- 


d 2 (pj{T*r, t*) 


dT* 2 


+ Rqy* 2 [F(t)e-^ Un + 4ir<p*JiT* v , t*)l 


(24c) 


If there is no visible source, the last term on the right hand side 
of (24a) or (24c) vanishes. Other conditions on the solutions of 
these equations are that < <p* v , <p* T , and T* must be finite and 


<P*iir* T , t) T*\t* t , t). 


(25) 


2.4. Physical Meaning and Magnitudes 
of Parameters 

The diffusion time may be written t D = pCA.TL Er /(khT/ 
L et ), where AT is an arbitrary, small rise in temperature and 
L ET = Ft' is the thermal extinction length. From this it is 
seen that t D is time to raise the temperature of a layer of 
thickness L ET by AT when a flux of heat kAT!L ET is con- 
ducted through the layer. 

The physical meaning of the parameter q = v a T*lTrE* T kT r 
may be seen as follows. The numerator, u^T^hx = J v , char- 
acterizes the radiant power incident in the medium, while the 
denominator, EpcT,. = kT r /L* ET , characterizes the flux of 
heat conducted by a temperature gradient of TJL* et . Now, 
the first term on the right-hand side of (24a) describes, in 
reduced form, the net flux of heat conducted through a layer 
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and the second term describes the net flux of power carried by 
thermal radiation through that layer. Hence irq is a measure of 
the ratio of the heat conducted by radiation to that carried by 
solid state conduction. 

To illustrate the magnitude of some of these parameters, 
ignore the subscripts v and T for the moment. Consider a 
simple medium composed of N particles per unit volume of 
equant, isotropically scattering particles of average diameter d 
and solid density p s . The bulk density is p — Np x Trd 3 /6, and 
the particle cross-sectional area is <r — ird 2 /4. Hence E = 
naQ E — (3/2 )4>(Q E /d), where c p = p/p s is the filling factor. 
For particles larger than the wavelength relatively close to each 
other in a regolith, Q E — 1 [Hapke, 1993a] so that the extinc- 
tion length is L E = FU 1 2d/34>. In loosely packed partic- 
ulate media, such as the upper surface of the lunar regolith, 4> 
can be as small as 0.3 [Carrier et al., 1973], so that L E ~ 2d. 
Apollo lunar soil samples have mean particle sizes around 50 
/xm [McKay et al., 1974]. Hence, for impact-generated rego- 
liths, L e ~ 100 /am, and E ~ 1 X 10 4 m~\ 

For a medium illuminated by the Sun, J v - 1360 ID 2 W 
m~ 2 , where D is the distance to the Sun in AU, so T r = 
(i tJ v /o- 0 ) II4 /VD = 524 /VD K. The temperature-indepen- 
dent component of the thermal conductivity of basaltic pow- 
ders, including lunar soil, depends on the density, but k — 
0.001 W m“' KT 1 is typical for densities around 1 g/cm 3 
[Fountain and West, 1970]. Hence q — 0.26 D 3/2 , which 
ranges from 1.07 for Mercury to 0.023 for a satellite of Jupiter. 

The specific heat per unit mass of lunar soil depends on the 
temperature, but averages about C — 550 J kg -1 K -1 over the 
temperature range from 100 to 350 K [Horai and Simmons , 
1972], Taking p - IX 10 3 kg m~ 3 , gives t D — 5.5 s. 


3. A Medium Heated From Below 

3.1. Solution 


The first solution that will be obtained is for a system con- 
sisting of a thick layer of powder in vacuum lying on a plate 
held at constant temperature T p and radiating into cold space 
from its upper surface. This might represent the situation in a 
laboratory measurement of emissivity where the surface of the 
sample sees only other surfaces at temperatures close to abso- 
lute zero. 

For this problem, dT*ldt* — 0, J „ — 0, <p* = 0, and q = 
a 0 T 3 / 7 rE* r k. Choose the reference temperature T r to be the 
temperature of the heated bottom plate T p . Let the emissivity 
of this plate be s p and the optical thickness of the powder at 
the plate be r Tp , so that its reduced optical thickness is r* Tp = 
Ct t t/>' Then equation (24c) is 


d 2 T*{ r* r ) d 2 <p* T ( r* T ) n 
dr* T 2 + 77(1 dr*, 2 


(26) 


where T* = T/T p is the temperature relative to that of the 
plate. 

Integrating (26) from 0 to t* gives 


dT* 

Jr% 


(t*) 


dT* 

~dr% 


(0) + p,q 


~d<p* T 


(t*t) 



= 0. 


Applying the boundary conditions (23c) and (24b) at the sur- 
face, this equation becomes 


dT* 


(t*) + Trq 


d<p* T 


(t* t ) = 2 77 q<p* Ts , 


(27) 


where the subscript s denotes the value at the surface: 

<P 7s = <P r(0)‘ 

Integrating again gives 

T*(t* t ) ~T*+ 7rq[<p%{r%) - <p * Ts ] = 2nq<p% s T%, (28) 


where 


T* - T*{ 0). 

This equation shows that the sum of the reduced temperature 
and 77 q times the reduced thermal flux increases linearly with 
optical depth through the powder. 

At t* t = T * Tp , T*( r* Tp ) = T* = 1. The boundary condition 
on <p* T comes from requiring that the upward radiance at the 
plate is the sum of radiance emitted by the plate plus the 
downward radiance from the powder reflected by the plate; 
that is, 

It i ( r r P ) = e/r 0 7; 4 /7r+ S p I T2 (j Tp ) , 

where S p is the reflection coefficient of the plate integrated 
over the upward hemisphere. By KirchhofFs law, S p = 1 - e p . 
Converting to reduced quantities, and using (14) and (23b), 
this boundary condition becomes 

<Pt( t *t p ) + A<p*( T* Tp ) = Sp + (1 - e p )[<p*( T* Tp ) - A<p*(t* p )], 
or 


57? ( A) = 2 ^; [i - A)]. (29) 

However, if the optical thickness of the powder is so large 
that near the plate the medium is in thermodynamic equilib- 
rium, then <p* T =* T* 4 = 1, so that at t* t = r* Tp , equation (28) 
becomes, 

1 ~ T*+ Trq{\ - <p*s) = 2nq(pr s T* Tp , 
which gives the simpler condition 


. 1 - T s + nq 

(PTS= 777/(1 +2r* rp )' 


(30) 


In order to complete the solution another relation between 
T* and (p* T is needed. This is provided by the solution to the 
thermal radiation equation (23a). Unfortunately, this equation 
is nonlinear and must be solved numerically. A program was 
written to solve (23) and (28) using the method of finite dif- 
ferences. The program starts by estimating a value for T 4 and 
using (30) to calculate a corresponding value for <p* Ts . Using 
the boundary conditions (23c) and (24b), the program inte- 
grates (23a) and (28) to find values of ip* T and T* for increasing 
values of t* t . If the guess for T* is incorrect, tp* T and T* 4 
rapidly diverge. Then a new value for T*. is chosen, and the 
process repeated until T* 4 and <p* T converge and also are 
insignificantly different from 1.00 at C T = r* Tp . 

Figure 2 shows the results of such a calculation for r* Tp = 
10, q = 0.26, and w* T = 0.36. This value for w * T is repre- 
sentative of the single scattering albedo of a particle larger 
than the wavelength in the thermal IR where strong rest- 
strahlen bands occur, and q is appropriate for an Earth labo- 
ratory simulating lunar surface conditions. The reduced tem- 
perature T* starts at T* = 0.6 with zero slope, as required by 
the boundary condition, and then rises with a steadily deceas- 
ing slope to 1 at i> /? . The reduced thermal flux increases 
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9 Ts ~ 


y t 


1 + 7: 


TV, 


(32) 


which is a second relation between (p Ts * and T* and can be 
used along with (30) to solve for these quantities. 

Using the boundary layer approximation, (28) becomes 


~ T*+ 77 q 


T* 4 ( t* t ) ~^e- 

y t 


2yrn 


9 Ts 


2 Tjq (p 'V T 7'. 


(33) 


At Tj- = t yy,, r* = 1, and this is 


1 + 7t * 

1/4 

* 9 Ts 

+ irq 

L 7 T J 

_ 


9 7 j 
7* 


,-2yrrr„ 


9 7j 


Figure 2. Reduced temperature and thermal radiance as a 
function of reduced optical depth in a layer of powder in 
vacuum lying on a heated plate and radiating into a cold 
environment from its upper surface. Parameters are t* T/ , = 10, 
w* T - 0.36, q = 0.26. Solid line is exact calculation; dashed 
line is boundary layer approximation. 


rapidly from its surface value q>* Ts — 0.07 with the required 
initial slope of 2 <p* Ts and then increases more gradually to 1 at 
T* Tp . When 2y* T T* r > 3, cp* r is within a few percent of T * 4 , 
and the medium is in approximate thermodynamic equilibrium 
with the thermal radiation. 

The rise in temperature near the surface has been pointed 
out previously by Henderson and Jakosky [1994], It is caused by 
the leakage of thermal radiation from the surface, which is also 
responsible for the near-surface gradient of q>* r . Within this 
boundary layer, which is of the order of l/y * T in optical thick- 
ness, the heat carried toward the surface by conduction is 
converted entirely into thermal radiation. 


3.2. Boundary Layer Approximation 

It would be convenient if an analytical solution to (23a) 
could be found. This is difficult because the equation is non- 
linear. However, an approximate analytical solution that is 
accurate to a few percent can be derived as follows. 

Rewriting (23a), 


= T* 4 


4 d 2 d% 

*4 TT * 

y t 


Now, the deep interior is in thermodynamic equilibrium so that 
<p* T — T* 4 there, showing that the second term on the right- 
hand side of the last equation is appreciable only near the 
surface of the medium. However, at the surface, dT*/dT* r — 
0, so that the temperature is almost constant there, T * = T*. 
In the near-surface region the solution to (23a) with T* con- 
stant and that satisfies the boundary condition (23c) is 


<P t(t 7) = TV 


<Pn 

* 

y t 


-2 yrrr 


These arguments suggest that the thermal flux may be approx- 
imated by 

<p*( T *) - - *1: e~ 2 ^. (31) 

y r 

At r = 0, this equation is 


n ^ * 

= 2 Trq ip Ts t T p, 


(34) 


which is a transcendental equation that may be rapidly solved 
for <p*. y by iteration. The solution is then completed by insert- 
ing cp* Tx into (32) and solving for T*. These approximate so- 
lutions (equations (31) and (33)) are plotted as the dashed 
lines in Figure 2. 

Converting (33) to real units gives the approximate solution 
for the temperature 


T = T, + 


q O 72 ^4 , .IT 

2<Pts£t t t —T+ <p Ts e 
0-0 . 77 7r 


—2{ryrrr 


+ <PTs 


(35) 


where T s ~ T r T* and <p Ts = (cr 0 /7 t)T 4 <p* T x - 

The flux of heat carried by radiative conductivity can be seen 
explicitly as follows. Substituting the boundary layer approxi- 
mation for cp* T into (27) for the gradient of T * gives 


dT* 

1 7, 


(t *) + Trq 


dT * 

4T*\tV) dr j (t*) + 2c p\e 


2y nr 


= 2TrqipV, 


Converting to unstarred quantities, this equation can be put 
into the form 


, dT(z) 


4 a 0 T 3 (z) d.T(z ) 
E t dz 


+ 2 tt£ 2 t < p T ,e ~ UrrnT 


= 2tt£ 2 t <Pt s . 


The first term on the left-hand side of the last equation is the 
flux of heat carried from the hot plate to the surface by con- 
ductivity, while the term in brackets describes the flux of heat 
carried by radiation. This expression shows that far from the 
surface, the radiation contribution to the thermal conductivity 
is proportional to T 3 , as is commonly assumed. However, this 
assumption is not valid in the boundary layer because it ne- 
glects the second term in brackets, which is important near the 
surface. The IR radiance that is measured by instruments 
above the surface comes from this boundary layer. 

The thermal flux radiated from the surface and the effective 
temperature at which the surface radiates into the space can be 
estimated using the boundary layer approximation and equa- 
tions (14) and (32). The thermal power emitted from unit area 
of the surface is 
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Pr — 


it/ 2 


I Tl (0) cos e dkl = / T ,(0) cos e 2tt sin ede 


= 77/7-1(0) = 7 t[<P 7 -( 0 ) + Acp r (0)] — 2 77 <p 7"(0) 

2y r 


~ ^ rp4 y T qr*4 _ 1 

= 277 - T r 9n - 2a 0 T r jT}—* T. ~ 1 + y 


<r Q Tt 


= e k * 0 n, ( 36 ) 

where s h is the hemispherical emissivity. Hence the medium 
radiates into space at an effective temperature T s , which is 
approximately equal to its actual surface temperature, and with 
a hemispherical emissivity 

£/j = 277/(1 + 7 *) = 2y t /(Ct + 7r)- ( 37 ) 


4. Planetary Regolith in Equilibrium 
With Sunlight 

4.1. Solution 

The second case to be considered is a semi-infinite particu- 
late medium permanently facing the Sun. In the classical 
steady state problem all absorption and emission of radiation is 
assumed to take place at the surface, so that the heat equation 
is simply d 2 T(z)/dz z = 0. The only solution to this problem 
is T(z) = const. If the temperature of the medium is T(z) = 
T e , and letting the medium have hemispherical albedoH,, and 
hemispherical emissivity s h , where both A h and e h are as- 
sumed to be constants of the surface independent of i and e, 
then the visible sunlight absorbed, (1 - A h )J v p 0 , is equal to 
the thermal power radiated, s h a 0 T 2, giving 

T c — [(1 — Ai,)J „/A(/£/i<ro ] ]l> ■ 

A more accurate model retains the assumption that the 
temperature throughout the medium is constant, but solves the 
visible and IR equations of radiative transfer to obtain A h and 
s h . This is done by Hapke [1993a, b], where it is shown that a 
medium of isotropic scatterers has a hemispherical albedo 
(directional-hemispherical reflectance) 

A h (w„ no) = 1 - 7 „H{w,„ p „), (38) 


giving 


E h (w T ) — 


2 7 r 
1 + 7 T 


I 


1 1 ~ 7 r \ 

6 1 + 77 / ' 


Thus, 


T e = 


J_JH JyH(W v , /X 0 ) 

cr o 2y t H { {w-,) 


1/4 


J v Pa 1 + 2pg 1 + 7r 1 

«j () 1 + 2y v p 0 2y r 1 1 ~ 7r 

6 1 + 7 r 


(44) 


(45) 


If the IR reflectance of the medium is small, as is often the case 
because of the strong reststrahlen bands there, y T is close to 
one and the factor (1 - 7 r )/6(l + 77-) in (44) and (45) may 
be neglected. 

The restrictions that the temperature of the medium is con- 
stant and that the particles scatter isotropically will now be 
dropped. Instead, the steady state heat equation and the visible 
and IR radiative transfer equations given in section 2 will be 
solved. 


For the steady state, dT*/dt* — 0, /x 0t ,«. - 
F(t*) = 1. Then (22)-(24) are 

1 d 2 (p* * w* v l + PJ 2p 0 * ^ 

4 dT7" J ' ><Pv 477 1 — p v e 

= const and 
++, (46a) 

1 d(p*, ^ w* p„ 
<Pw(0) “ 2 ~dr* ” 477 1 - p v ’ 

(46b) 

i Cl ip T .,* 21 ..* 

Cdv?- yA ‘ fr - T ’’ 

(47a) 

+<°) = 

(47b) 


3, (48a) 

dT * 

^ ( o, = o. 

(48b) 


and a hemispherical emissivity 

s h (w T ) = 2y T H l (w T ), (39) 

where H(w, x) is the solution of the integral equation 

w f'H ( w , x') 

H(w, x) = 1 + —xH(w, x) ~ x + x dx' , (40) 

7 0 


and may be approximated by 

H(w,x) - (1 + 2x)/{\ + 27*), (41) 

so that 

A„{w xo ju, 0 ) “ (1 - 7 .;)/(! + 27#o)> (42) 

and H x (w) is the first moment of the H function and may be 
approximated by 


Equation (46a) for the visible flux is independent of the 
other two equations and can be solved separately. The solution 
satisfying the boundary condition that the flux is finite at in- 
finity is 

cp*, = A'e + B'e (49a) 


The coefficient B ' can be found by substituting this solution 
into the differential equation (46a) and equating coefficients of 
exp r*J p* Ul ); A' can be found by substituting the solution 
into the boundary condition (46b). This gives 


A ' 


wl 4^(1 + PJ 2pt) 
47ra- ^(1-47^’ 


(49b) 


, __ 1 2 ixl,(l 

B 477 


7t)[(l + 2/xtJ + m + 2y* 2 p*J] 
(l-47>? B 2 )(l-i^) 7 




H(w, x) xdx 


1 + 7 


<«> 


Inserting this solution into (48a), the heat transfer equation 
becomes 
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However, T* must retain finite as r*r -> The only way this 
can be true is if the sum of the coefficients of t* t vanishes; thus 


2ir<p* Ts =y* v 2 (An* Qv +B/2y*). 


(55) 


After substituting (50b) and (50c) and simplifying algebra- 
ically, this becomes 

* * 

(pTs= 2ir(l + 2y*fi*j 

(1 + 2fiov) + PJ( 1 - 2yt- 2 7 :VoJ 

(1 - |3J 

Thus the solution to the heat transfer equation is 
T*(t*) - T*+ 7rq(cp*r( t*) ~ <PrJ 


(56) 


Figure 3. Reduced temperature and thermal radiance as a 
function of reduced optical depth in a regolith in vacuum in 
equilibrium with sunlight, showing the effect of varying the 
angle of incidence i. Parameters are w* = 0.36, w* T = 0.36, 
q = 0.26. Solid line is exact calculation; dashed line is bound- 
ary layer approximation. 


d 2 cp* r 


d 2 T* 

dTj 

where 


-Rqy* 2 [Ae~ RTf/lMt + Be 2y ' :Rjr ], (50a) 


A — 1 + 4 ttA 


1 + 2/itf (l-2gt-|3 s (l-2y:VL) 

1 - 4y?n% 1 - j8„ 


(50b) 


„ „ 2 J aL(l - 7t)[(l + 2^„) + J3„(l + 2 7 :Vo,)] 

(i - 4 T :vj)(i - pj 

(50c) 

Integrating (50a) once with respect to t* t from 0 to and 
applying the boundary conditions gives 




-Rrrljxa 


') + 5/4 (1 - e^ 2y ' r,iT )l (57) 


Now, deep in the interior of the medium as t*- — > <p* —> 

0, T* and (p* r become constant; hence d 2 q>* T ld r r\ = 0, and 
cp^(co) = t* 4 ( oo). Thus, as t — » co, equation (57) becomes 


T*(«>) ~T*+ t TqT*\°o) 

- q[ tt<p% s + (MR) (ytVo M + 5/4)]. 


(58) 


dT * , 

Ti% ^ + 


thpj , 

-h < T *) ~ 


*2 


B 

2y* v 


where 

<P* n = vX o). 

Integrating again gives 

T*(t* t ) ~ T*+ 'Trg[<p%(T* T ) - (pr x - 2(p* Ts T* T ] 

^ (e R ^ /lxl - 1) + T* 

K 


- -qy* 2 \Ag 


B 

2v* 

^ 1 V 


2y 


*R 


where 


( e -2y,tof_ 1) + T * 


r* = r*( 0 ). 


This is an alternate boundary condition. 

Equations (47) and (48) were solved numerically by the 
method of finite differences in a manner similar to that used in 
section 3. The procedure starts by guessing a value for T* and 
integrating inward from the surface, with <p* Tx given by (56). If 
the guess for T * is incorrect, the numerical solutions for (p* T 
and T* 4 rapidly diverge, in which case the integration is 
stopped, a new value for T* is estimated, and the procedure is 
repeated until a value is found for T*. that allows c p* r and T* 4 
to converge with increasing t* t . 

The results of the numerical calculation are shown as the 
solid lines in Figures 3-6, which illustrate the effects of varying 
the parameters of the solutions. In Figures 3-6 it was assumed 
that the particles are isotropic scatterers and that E v - E r , so 



(53) 


(54) 


Figure 4. Reduced temperature and thermal radiance as a 
function of reduced optical depth in a regolith in vacuum in 
equilibrium with sunlight, showing the effect of varying the 
visible single scattering albedo w*. Parameters are w* T = 
0.36, q = 0.26 , / = 0. Solid line is exact calculation; dashed 
line is boundary layer approximation. 
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that R = 1 and t* t = r*. The value for the visual single 
scattering albedo of w* = 0.36 might be representative of a 
relatively low albedo body, such as the Moon or an asteroid. 
The high value, w*, = 0.91, might be representative of some- 
what dirty ice on the surface of a rocky-icy satellite, and w* = 
0.9996, representative of clean, fresh frost or snow. In the 
thermal IR, w* T = 0.36, is probably representative of most 
materials. A value of w * T = 0.91 probably is unrealistically 
high for most materials of planetary interest because rest- 
strahlen bands are common in most materials at IR wave- 
lengths; however, it is included to show the effect of a material 
of high IR albedo, such as an alkali halide. 

In all cases, as the optical depth r* T increases, <p*r rises 
rapidly from a value of <p at the surface, as required by the 
boundary condition, and then levels off either to a constant 
value or to a more gentle rate of increase. The rapid rise is a 
direct result of leakage of thermal radiation from the surface. 
The medium approaches thermodynamic equilibrium when 
t* t > 3. 

The behavior of T* is more complex. The slope of T* is zero 
at the surface, as required by the boundary condition, but as t* t 
increases, T* may exhibit a positive or negative slope, depend- 
ing on the relative value of y* T compared with y* and (Xq. If 
these quantities are comparable, T* has a positive slope. How- 
ever, the slope of T* may be small or even negative if y* T is 
small. Heating by unscattered visible radiation occurs over an 
optical depth of the order of l/fio v = 1/^p-o an< l by multiply 
scattered visible radiation over 1/y* = Z,Jy v , but the IR flux 
is radiated from the surface over an optical depth of the order 
of 1/y r = Ct/Jt' If 1/Y r is larger than l/fi* 0v , or l/y*„ then 
heat must be conducted from the visible heating source, which 
is closer to the surface, to supply power to the deeper IR 
radiative sink. This requires a negative temperature gradient. 
Thus, when i is large a negative greenhouse effect can occur. 
Henderson and Jakosky [1994, also submitted manuscript, 1995] 
have previously pointed out that large positive subsurface tem- 
perature gradients are possible, but they did not investigate 
cases where the gradient was small or negative. It should also 
be emphasized that the widths of the gradients in actual depth 



Figure 5. Reduced temperature and thermal radiance as a 
function of reduced optical depth in a regolith in vacuum in 
equilibrium with sunlight, showing the effect of varying the 
thermal single scattering albedo w* T . Parameters are w* = 
0.36, q = 0.26, i = 0. Solid line is exact calculation; dashed 
line is boundary layer approximation. 



Figure 6. Reduced temperature and thermal radiance as a 
function of reduced optical depth in a regolith in vacuum in 
equilibrium with sunlight, showing the effect of varying* the 
radiative conductivity parameter q. Parameters are w* = 
0.36, w* T = 0.36, i = 0. Solid line is exact calculation; 
dashed line is boundary layer approximation. 


z (rather than in reduced depth T* r ) will depend inversely on 
the particle density or filling factor. 

As shown in Figure 6, the parameter q has a relatively small 
effect on the reduced variables (although, of course, the abso- 
lute effect can be large). The reason can be seen from (58) for 
the total rise in temperature. As q increases, the visible heating 
increases, which increases r*(°°) — T *, but the radiative 
conductivity also increases, which limits the rise in tempera- 
ture. 

Substituting the boundary layer approximation 

* 

<pUtt) ~ T*\r%) - c~ 2rM , (59) 

y t 

into (57) gives an approximate expression for T*, 

= T*~ TrqT*\T* T ) - X <p% s e~ 2yU * - <p* Ts ] 
y t 




ytVk 4(1 -e 


a -R tt/P-U, 


B 


) + Td-* 


_ 0 -2y„RTT\ 


(60) 


where 


r * = 

S 



1/4 


(61) 


and (p* Ts is given by (56). The boundary layer approximation 
(equations (59)-(61) is plotted as dashed lines in Figures 3-6. 
The exact and approximate solutions are within a few percent 
of one another everywhere. 


4.2. Emitted Thermal Radiance and 
Hemispherical Albedo 

The thermal flux radiated from the surface, the effective 
temperature at which the surface radiates into space, and the 
hemispherical emissivity can be determined in the same way as 
was done in section 3 for the case of the powder heated from 
below. From (36), the thermal power emitted from the surface 
is 
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Figure 7. Thermal flux emerging from the surface of a rego- 
lith in vacuum in equilibrium with sunlight as a function of the 
angle of emergence e for various angles of incidence i . Param- 
eters are w v = 0.36, w T = 0.36, j3„ = 0, p T = 0, q = 0.26. 
Solid line is exact calculation; dashed line is constant subsur- 
face temperature. 


Pr = 


7:n(0) cos e dCl = 2 tt(p t (0) — 2 7 rJ v (p* Ts 


(T 0 y t 
= 2 tt — 7 r ^— 

7T 1 + y T 


_ 




l + y 


* vj*. 


(62) 


Hence, to the accuracy of the boundary layer approximation, 
the effective radiating temperature of the surface is 


T s = T r T* = 


77 T 1 + yr * 

J v * <Pts 

cr 0 7 r 


and the hemispherical emissivity is 

s h = j(l + 7 j) — 2y T !{l T + 7 T ), 


(63) 


(64) 


which is the same result obtained in section 3.2 for the powder 
heated by a bottom plate. 

The visible hemispherical albedo (directional-hemispherical 
reflectance) Af,(p 0 ) of a surface in radiative equilibrium with 
sunlight can he found from (56) for <p* Ts by noting that the 
thermal power emitted per unit area must equal the visible 
power absorbed per unit area; that is, 


Pt~ 2ttJ v (p* Ts — / r ,p, 0 [l A h (po)]. 


Thus 

A h (p<)) = 1 


2 tt 


* 1 

VTs = 1 


vl * 


Pu • " (1 + 27VoJ 

(1 + 2p*J + fi v (l - 2y* v - 2y**pt) 


1 - Pv 


- 1 - 


SvVv 


1 + 2 £„y „p 0 

(1 + 2 ^ 0 ) + J3„(l - 27 , U v - 2ylfi 0 ) 

1 - P v 

which can be further simplified to 


(65a) 


L~y v L~p v y v 

1 + 2f„7 v p 0 1 - p v 


(65b) 


Comparing (65 a) and (38) for A h shows that the two are 
identical when p v = 0. Comparing (64) and (39) for s h shows 
that when p T — 0, the two are the same if H } (w T ) is replaced 
by its low-albedo approximation H^(w T ) - 1/(1 + 7 T ). This 
approximation is probably valid in most cases of interest, be- 
cause strong reststrahlen bands cause the IR albedos of most 
substances to be small. 

By symmetry [Hapke, 1993a], the visible hemispherical- 
directional reflectance has exactly the same form as (65a) or 
(65b), except that p 0 is replaced by p. Also, the IR directional- 
hemispherical, and hemispherical-directional reflectances have 
the same forms as their visible counterparts, except that sub- 
scripts v are replaced by T. 

Hapke [1993a] shows that the IR radiance emerging from 
the surface of a particulate medium is given by 


It(p 05 p) — 


w T 

Att 


I T ( t t , p, 0,')p T (Q', H) dfl' 


+ y 2 r ~T 4 (r T ) 

7 T 

In the two-stream approximation with the hemispherical asym- 
metry factor, the integral over solid angle is 


T7-//A 


d t t 


( 66 ) 


w T 

Att 


I t( +7 1 , pj f 1 )py(fl , fl) d£l — Wj(Pf{ Ty) 


+ Pfw . (6 7) 

Substituting this into ( 66 ), changing to reduced quantities, and 
using (23b) gives 


I rip o> p) — J v 


1 - 


-„V r ( T* r ) + 21 


1 & 


I3 t W : 


d<p%{T%) 

dr* T 


+ ttyVVir * T ) 


* ■ 
ll T 


( 68 ) 


Equation ( 68 ) was evaluated numerically to give I* T (p {) , p) 
= I T (p 0 , p)U v versus e for a lunar-like regolith. The results 
are shown as the solid lines in Figure 7. 


4.3. Directional Emissivity and the Constant Subsurface 
Temperature Approximation 

The integrand in ( 68 ) contains terms proportional to <p* T , 
d(p* T /dr* T and T * 4 multiplied by exp (-t * T /p* T ). This inte- 
grand is appreciable only over a distance of the order of 1 /p* T 
from the surface. However, because of the zero initial slope of 
T*, within the layer, T* is nearly constant. Hence the boundary 
layer approximation with T* = T * may be used to obtain an 
approximate evaluation of ( 68 ). This gives 


I T (pa, p) - 


= jjt 


lry*r 

IT 2 7 VI- 


(1 + 2/4) + /3 r (l - 27* - 27 *7+4) 
1 - p T 


(69) 


Writing this in the form 
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n) = ^ r; (70) 


shows that the effective radiating temperature of the medium 
is approximately 



It + 7t 
It 



1/4 


(71a) 


where <p* Ts is given by (56) and the directional emissivity is 
given by the approximate expression 

£tJt (1 + 2£^ju,) + /M 1 — 2y T /£ T ~ 2y r^) 

6rf(/A) “ 1 + 2£ r y£ 1 - ^ ' 

(71b) 


If the particles scatter isotropically in the IR (fi T - 0), the 
emissivity is 

sAv) = 7 t yrH(w r , /x), (72) 

which has been previously derived by Hapke [1993a]. 

It should be emphasized that (56) for <p* Ts and (65a) or (65b) 
for A h (^o) and the visible and thermal-IR directional- 
hemispherical and hemispherical-directional reflectances are 
exact and do not depend on the validity of the boundary layer 
approximation. Equations (64) for s h and (71a) for T s depend 
on the boundary layer approximation but not the constant 
subsurface temperature approximation. Expression (71b) for 
s d (p) depends on the validity of both the boundary layer and 
constant subsurface temperature approximation. All of these 
expressions include the effects of the gradient in cp* T , which is 
closer to the surface than the gradient in T* . 

By combining (56), (71a), (71b), and (72), the surface tem- 
perature can be written in the form T s - [(1 - A h )J v p 0 l 
e/i o- 0 ] 1/4 , which is identical to the expression for the radiative 
equilibrium temperature T e derived at the beginning of this 
section under the simple assumption that the temperature of 
the medium is constant. Thus, allowing for a solid-state green- 
house increases or decreases the temperature distribution in 
the interior of the medium but does not alter the surface 
temperature. 

The approximate expressions (69)-(71) for/ r (p. 0 > P-) versus 
e were calculated for the same values of the parameters as used 
in the numerical evaluation of (68). The results are shown as 
the dashed lines in Figure 7. The agreement between the exact 
and approximate values is good when either i or e is large, but 
there is a discrepancy of as much as 5% when both i and e are 
small because of the neglect of the subsurface temperature 
gradient in evaluating the integral in I T (p o, M-)- H e is ^ ar S e > 
the integrand in the expression for I T is appreciable only close 
to the surface, where the temperature is approximately con- 
stant. If i is large, the temperature gradient is small and its 
neglect introduces little error. The gradient influences the 
emergent radiance only when the surface is both viewed and 
illuminated at nearly vertical angles. 

The hemispherical-directional reflectance in the thermal-IR 
is given by (65) with replaced by p and subscript v replaced 
by T. By Kirchhoffs law [Hapke, 1993a], the directional emis- 
sivity is the complement of the hemispherical-directional re- 
flectance. Comparing (71b) with (65a) shows that this is indeed 
true. However, (71b) is an approximation which is valid only 
when a large subsurface temperature gradient is not present. 
This is in agreement with the result previously obtained exper- 



Figure 8. Directional emissivity s d (p) as a function of angle 
of emergence e and the complement of the hemispherical 
albedo 1 - A h { ju. 0 ) as a function of angle of incidence i for 
different values of the thermal or visible single scattering al- 
bedo w and hemispherical asymmetry parameter (3. 


imentally by Salisbury et al. [1994], who showed that KirchholFs 
law is violated in media with large near-surface temperature 
gradients. 

The effect of single scattering albedo and asymmetry factor 
on e d (p) is plotted for a high and a low albedo in Figure 8. 
Because 1 - A h ( p 0 ) has the same functional dependence on 
I as e d (p) does on p, these curves also describe 1 - 
A h (i x 0 ), as indicated in Figure 8. Making the particles more 
forward scattering increases the directional emissivity and de- 
creases the hemispherical albedo but does not change the 
shapes of the curves appreciably. 

4.4. Spectral Boundary Layer Approximation 

For future use, it is of interest to write down explicitly the 
unreduced boundary layer approximation for the thermal and 
spectral fluxes in the IR. Converting the boundary layer ap- 
proximation (equation (59)) to unstarred quantities gives for 
the radiative flux integrated over thermal wavelengths, 

vAt r) “ - n t t ) - £ <p^- 1Crn "- (73) 

IT Jt 

By analogy, the boundary layer approximation for the radiance 
averaged over angle in the thermal-IR wavelength region 
where the visible source spectral irradiance is negligible is 


where 


f.K) = l VAT) - ^ 

TT Ta 


(74) 


<P\S 


7a 1 

-U k (T s ), 


L + 7a tt 
T s = T r T 


(75) 


T* is given by (61) and (56), and T r is the blackbody radiative 
equilibrium temperature. 


5. Conclusions 

The temperature and thermal radiance in a planetary rego- 
lith are governed by three equations: the equation of radiative 
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transfer for visible and for IR radiation and the heat transfer 
equation. These equations were written in terms of the two- 
stream approximation, including anisotropic scattering, which 
allows numerical solutions to be rapidly obtained. An approx- 
imate analytic solution, the boundary layer approximation, to 
the nonlinear thermal radiative transfer equation was ob- 
tained. This approximation allows an analytic solution to the 
heat equation to be derived that gives the distribution of tem- 
perature and thermal radiance as a function of depth, and that 
is accurate to a few percent. The equations intrinsically contain 
radiative conductivity and result in a changing temperature 
with depth, i.e., a solid-state greenhouse effect. 

Using this formalism, the steady state cases of a powder 
heated form below in the laboratory and a planetary regolith in 
equilibrium with sunlight were studied. Subsurface tempera- 
ture gradients that are positive, negative, or zero can occur, 
depending on the thermal and radiative constants and the 
geometry. Large gradients may occur in a layer only a few 
particles thick or may be spread out over a large depth, de- 
pending on the visible and IR radiative constants of the me- 
dium. Approximate analytic expressions for the directional- 
hemispherical albedo (equation (65)), the hemispherical 
emissivity (equation (64)), and the directional emissivity 
(equation (71)) for regoliths of anisotropically scattering par- 
ticles were derived. This model was also used to study three 
interesting problems in the thermal emission from planetary 
regoliths by Hapke [this issue]. Using the same formalism, an 
improved expression for the bidirectional reflectance of a re- 
golith of highly anisotropically scattering particles was also 
derived and will be the subject of a separate paper (B. Hapke, 
manuscript in preparation, 1996). 

Notation 

A fl hemispherical albedo, equal to directional- 
hemispherical reflectance. 

C specific heat per unit mass. 
c 0 speed of light in vacuum. 
dj mean diameter of a particle of type /', equal to 
(4oy/ir) 1/2 . 

E x (z), £„(z), E t (z ) spectral visible, and thermal volume 
extinction coefficients, equal to 
l j N j (z)c r jQ EKj , 2 j N j (z)a j Q EvJ , and 
2jNj(z)<TjQ EvTi respectively. 
e zenith angle at which radiance emerges from 
the medium toward the detector. 

% X ,% T spectral and thermal volume emission 

coefficients, equal to K x (z)/E x (z) — y 2 , 
K t (z)/E t (z) = y\ 

F(t) function describing the time dependence of 
the irradiance. (For a rotating body F(t) is a 
square wave with period P: F(t) — 1 when 
0 < t < P/2, and F{t) = 0 when P/2 < t < P.) 
G x (z, XT, XI), G v (z, CL', XI), G t {z, CL', CL) spectral 

visible, and thermal volume angular scattering 
coefficients equal to Xy /V ; (z) o'jQs^p^Cl' , XI), 

2 j Nj(z)ar j Q Sv jp v j(CL ', XI), and 
2jNj(z)o-jQ STj p rj (CL' , XI) respectively, for 
scattering radiation traveling in direction XI' 
into direction XI. 

g phase angle, equal to angle between the 
directions to the source and detector as seen 
from the medium. 


g' angle between the direction tt — XI' from 

which the radiance comes and the direction XI 
into which it is sca ttered. 

H(w, x ) (1 + 2x)/(l + VT~ w x). 

H x {w) first moment of the H functions. 
h 0 Planck’s constant. 

I x (z, ft, t), I v (z, ft, t), I T (z, ft, t) spectral, visible, and 
thermal radiances (diffuse radiant power per 
unit area per unit solid angle), respectively. 

I vl (z, t ), I T \(z, t ) hemispherically averaged visible and 
thermal radiances traveling into the upward 
hemisphere. 

I v2 {z, t) I T2 (z, t ) hemispherically averaged visible and 

thermal radiances traveling into the downward 
hemisphere. 

i(t) angle between ft 0 and z. 

J x , J v spectral and visible irradiance, respectively 
(collimated radiant power per unit area 
perpendicular to the direction of propagation) 
incident on the upper surface of the medium; 
for planetary applications, J x has the spectrum 
of a blackbody with a temperature of 5770 K, 
and J v = if ID 2 , where if is the solar 
constant, if - 1360 W m“ 2 and D is the 
distance to the sun in astronomical units. 
j subscript denoting particle type. 

K x (z), K v (z), K T (z) spectral, visible, and thermal volume 
absorption coefficients, equal to 
l J N J (z)a j Q Akj , ljNj(z )<JjQ Avj , and 
l j N j (z)a j Q AT . respectively. 
k solid state thermal conductivity. 
k 0 Boltzmann’s constant. 

F Exi , L et visible and thermal extinction lengths, equal to 
1/E V and 1/E r respectively. 

Nj(z ) number of particles per unit volume of type j. 
P period of rotation of the body. 

p Xj (CL' , ft), p vj (CL', ft), g r/ (fl', ft) spectral, visible, and 
thermal angular phase functions of a particle 
of type j, respectively; if the particles are 
randomly oriented, the particle phase 
functions depend only on g ' . 

p x {z, ft', ft), p v (z, ft', ft), p T (z, ft', ft) spectral, 
visible, and thermal average volume particle 
scattering functions, equal to G x (z, ft', CL)/S x (z), 
G v (z, ft', ft)/S„(z), and G T (z, ft', CL)/S T (z). 

QaxA z )> Qavj(z), QatA z) spectral, visible, and thermal 
volume absorption efficiencies of a particle of 
type /. 

QexA z )> Qe v A z )’ QetA z ) spectral, visible, and thermal 
volume extinction efficiencies of a particle of 
type/. 

QsxA z )’ QsvA 2 )’ QstA z ) spectral, visible, and thermal 
volume scattering efficiencies of a particle of 

type/. 

S x (z), S v (z), S T (z ) spectral, visible, and thermal volume 
scattering coefficients equal to HjNj(z)ajQ sxj , 
'F j N j {z)(j j Q Svj , and 2 y /V ; (z)oy0 iS7y 
respectively. 

T(z, t) absolute temperature. 

T s (t) surface temperature, equal to T(0, t) 
t time. 

U X (T) Planck function. 
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w A , w v , w T spectral, visible, and thermal volume-average 
single scattering albedos, equal to 
S \(z)/E x (z), S v (z)/E v (z), and 
S t (z)/E t (z). 

z vertical position (positive upward). 

/3 A , 0„, Pt spectral, visible, and thermal hemispherical 

asymmetry scattering parameters, respectively. 
Tt spectral, visible, and thermal albedo factors 
equal to, (1 - tv A ) 1/2 , (1 - w t ,) l/2 , and 
(1 - w T ) 112 respectively. 

<p A (z, t), <p v (z, 0’ 0 spectral, visible, and thermal 

average isotropic radiant fluxes, respectively. 

<P Xs (t), <p«*(0> <pts ( 0 <Pa(o, 0. <p w (0» 0> and ‘Pt-CO, 0 
values of the spectral, visible, and thermal 
fluxes equal to, respectively, at the surface. 

Ac p A (z, /), A<p„(z, t), A<p r (z, t) spectral, visible, and 
thermal average radiant up-down flux 
differences, respectively. 

<fi filling factors, equal to p/p s . 
ft angle between LL and the normal to the 
surface. 

A wavelength, 
p, cos e. 
ju 0 (t) cos i(t). 
p bulk density. 
p s particle solid density. 
cTy(z) mean geometrical cross-sectional area of a 
particle of type j. 
cr () Stefan-Boltzmann constant. 

t a (z)=J7£a( 2 ') dz > T r’( z ) = JT E X Z ') dz ’ 

t r (z) = £ r (z') dz spectral, visible, and 

thermal optical depths, respectively. 
dr k -£ a (z) dz, dr v = ~E v {z) dz, 
d T r = -S^z) dz. 

fl direction into which the radiance is traveling 
through an element of volume within the 
medium. 

a 0 direction toward the source from the medium 
(the incident irradiance travels into the 
direction i r — Il 0 ). 

£a=(1 - P\ w xY /2 ’ L = (1 - 0»”O ,/2 > 

= (1 - j3 r w r ) 1/2 , spectral, visible, and 
thermal asymmetry factors, respectively. 

The following relations also hold: 

j ~~ Qa\ j e Tj QaTJ- 

«a(z) = K x (z) % t (z)=K t {z). 

Qs\j + G/l A, - QEKp Qsvj + 0 /lrv = QEvp QsTj + 0 / 17 ) = QeTJ- 
S, + K k = E x , S v + K v — E m S t +K t =E t . 


yl = KJE X , y 2 v - KJE V , y\ - K T IE T . 

If the filling factor 4> of the medium is not small, Nj(z) must be 
multiplied by the factor -ln(l - (/>)/<#>. 
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